Barrier crossing driven by Levy noise: Universality and the Role of Noise Intensity 
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We study the barrier crossing of a particle driven by white symmetric Levy noise of index a and 
intensity D for three different generic types of potentials: (a) a bistable potential; (b) a metastable 
potential; and (c) a truncated harmonic potential. For the low noise intensity regime we recover the 
previously proposed algebraic dependence on D of the characteristic escape time, T csc ~ C{a) /D^ a ' , 
where C(a) is a coefficient. It is shown that the exponent n(a) remains approximately constant, 
H ~ 1 for 0<Q<2;ata = 2 the power-law form of T esc changes into the known exponential 
dependence on 1/D; it exhibits a divergence- like behavior as a approaches 2. In this regime we 
observe a monotonous increase of the escape time T csc with increasing a (keeping the noise intensity 
D constant). The probability density of the escape time decays exponentially. In addition, for 
low noise intensities the escape times correspond to barrier crossing by multiple Levy steps. For 
high noise intensities, the escape time curves collapse for all values of a. At intermediate noise 
intensities, the escape time exhibits non-monotonic dependence on the index a, while still retaining 
the exponential form of the escape time density. 

PACS numbers: 05.40.Fb,02.50.-Ey,82.20.-w 



I. INTRODUCTION 

In his seminal paper [1], Kramers proposed to model 
chemical reaction rates as the diffusion of a Brownian 
particle, initially located in a potential well, across a po- 
tential barrier of finite height. Meanwhile, Kramers' the- 
ory has been applied to a much more general range of pro- 
cesses associated with the barrier crossing of a physical 
entity experiencing random kicks fuelled by its contact 
to a thermal bath [2-5] . 

Since Kramers' solution numerous different ways have 
been reported to access the barrier crossing problem, in 
particular, to find the mathematically most convenient 
formulation, compare, for instance, Refs. [3-6]. The re- 
sult for the characteristic escape time reads [7] 
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The barrier crossing in the presence of long-range mem- 
ory effects was, inter alia, modeled in terms of a gener- 
alized Langevin equation [12, 13], or via a subdiffusive 
fractional kinetic approach with Mittag-Leffler survival 
[14, 15]. Subdiffusion is usually associated with a long- 
tailed waiting time distribution ip(t) ~ A 7 r 7 /i 1+7 with 
< 7 < 1 rendering the resulting continuous time ran- 
dom walk process semi-Markovian, its hallmark being the 
power-law time-dependence (x 2 ) cx Kjf of the mean- 
squared displacement in absence of an external potential 
[8, 9, 16]. 

While in subdiffusion the waiting time between succes- 
sive jump events becomes modified such that the mean 
waiting time J Q tijj(t)dt diverges and consequently no 
natural time scale separating microscopic and macro- 
scopic events exists, the distribution of the lengths of in- 
dividual jumps is narrow. The converse is true for Levy 
flights: Here, the lengths of the jumps are distributed 
according to the long-tailed jump length distribution 



where U (x) is the dimensionless potential, a prime stands 
for the derivative with respect to the coordinate x; x m i n 
and ir max are the points, where the potential U (x) at- 
tains its minimum and maximum, respectively; and D is 
the diffusion coefficient (noise intensity) of the diffusing 
particle, stemming from its coupling to the heat bath. 
Eq. (1) is based on the assumption that the barrier is 
high (or cquivalently, the noise intensity is low) , and there 
exists a constant probability flux across the barrier max- 
imum. 

Random processes in complex systems frequently vio- 
late the rules of Brownian motion. Thus, the presence 
of static or dynamic disorder might give rise to mem- 
ory effects causing subdiffusion [8, 9], and possible de- 
viations from standard exponential relaxation [10, 11]. 
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with < a < 2 [8, 15, 17, 18]. Thus, the variance 
x 2 X(x)dx of the jump lengths diverges. Such power- 
law forms of the jump lengths have been recognized in a 
wide number of fields [9]. Prominent examples for such 
genuine Levy flights are known from noise patterns in 
plasma devices [19], and from random walks of particles 
or excitations along a fastly folding polymer, where the 
walker is allowed to cross the small gap between two seg- 
ments of the chain, that are close by in the embedding 
space due to polymer looping [20] . In the latter case, the 
exponent a is in fact related to the critical exponents of 
the polymer chain. Further examples come from fluctua- 
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tions in energy space in small systems [21, 22], and from 
palcoclimatic time series [23]. 

From a mathematical point of view the occurrence 
of long-tailed distributions appears quite natural due 
to the Levy-Gnedenko generalized central limit theorem 
[24, 25]. Indeed, the tails of probability densities of the 
type of X(x), Eq. (2), are obtained from the characteris- 
tic function of a symmetric a-stable distribution of the 
form 



A(fc) = ,^{A(x)} 



f 

J —i 



e lkx \{x)dx = exp (-c|fc| Q ), 



( 3 ) 

where c > 0. The power-law asymptotics at large |x|, 
Eq. (2), appear immediately from the expansion of the 
characteristic function Eq. (3) in the limit of small k. 

An important question arises when replacing the Brow- 
nian particle in a barrier crossing process by a particle 
executing Levy flights. This situation can be modelled 
by a particle subject to Levy stable noise, on the level 
of the Langevin equation. First steps in this direction of 
addressing have been taken, as reported in Refs. [26-29]. 
In the present work, we report extended simulations re- 
sults for the usually studied case of low noise strength D, 
observing a pronounced step-like behavior of the depen- 
dence of the exponent /z on the Levy index a. We also 
explore the case of intermediate and high noise strength, 
finding a quite rich behavior in the parameter space, in- 
cluding an optimum a for the escape time. 



II. UNDERLYING LANGEVIN EQUATION 

We start from the Langevin equation for a particle em- 
bedded in an external potential field and subjected to a 
random noise, 



dx(t) 
dt 
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(4) 



where x(t) is the dynamic variable (particle position), m 
the mass, 7 the friction constant, and U(x) an external 
potential. The noise £ Q (t) is a white, symmetric a-stable 
noise. Eq. (4) is understood in the following way [30]. 
Integrating Eq. (4) over the interval [t, t + At] , we obtain 



x(t + At)-x(t) 



^ rt+At 

Jt 



dU(x(t')) 
dx 



where 



dt' + L atD (At), 
(5) 
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is an a-stable process with stationary independent incre- 
ments and characteristic function 



p L (k, At) = cxp (-D\k\ a At) . 



(7) 



D is the intensity of the Levy noise. With the use of 
Eqs. (6) and (7) it is straightforward to show [30] that the 



discrete time representation of Eq. (5) at times t n — nSt, 
n = 0, 1, 2, . . . for sufficiently small time step St is 



1 dU(x n ) 
7717 dx 



St+(DSt) 1/a Ui(n), (8) 



where {£0,1(77)} is a set of random variables possess- 
ing Levy stable distribution A(x) with the characteristic 
function (3) and c = 1. 



III. SIMULATIONS RESULTS 

Before addressing the barrier crossing problem in the 
presence of Levy stable noise analytically below, we 
present results from extensive simulations. In these sim- 
ulations, we employ three types of potential profiles: 
bistable, metastable, and truncated harmonic potentials 
(see Fig.l) defined as follows: 



(9a) 
(9b) 

(9c) 



Let us turn to dimcnsionlcss variables. To this end we 
substitute x — > x/x, t — > t/t, and D — ► D/D into the 
discrete time Langevin equation (8), and choose the ap- 
propriate constants for rcscaling, x, t, and D, for each 
potential type. Taking into account that £ Q ip/'t) = 
?" 1/a ^a(0 [30], we find 
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respectively. In terms of these rescaled variables the dis- 
crete time Langevin equation (8) reads 



dU{x n ) 
dx 



5t+(DSt) 1/a S; aA (n) 7 (11) 



where the potential U(x n ) refers to, respectively, Ui, Ui, 
or t/3 from Eqs. (9), with g 1 = g 2 = L = 1. 

In the simulations, for each type of the potential the 
particle starts from the bottom of the potential well. The 
adsorbing boundary is located: (i) for the bistable poten- 
tial in the saddle point x = 0; (ii) for the metastable po- 
tential far to the right of the saddle point, at x = 10; and 
(iii) for the truncated harmonic potential at the bound- 
aries, |x| = 1. The process defined by Eq. (11) is repeated 
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Figure 1: The three types of potential profiles Ui (bistable potential), Ui (metastable potential), and U3 (truncated harmonic 
potential) considered in the simulations; for gi, g%, L > 0. 



until the particle reaches, or crosses the adsorbing bound- 
ary; then the process in Eq. (11) is restarted. This proce- 
dure is performed 100,000 times, for fixed D and a, and 
the mean escape time is then calculated [31]. 

The mean escape times T csc of the diffusing particle 
as function of noise intensity D for the three potential 
profiles are shown in Fig. 2 for different values of a rang- 
ing from 0.1 to 2 [32]. It is clear that the curves for 
a < 2 obey a different law in comparison to their Gaus- 
sian counterpart. We observe three different regimes for 
the dependence of T csc on the Levy noise parameters, 
which can be classified by the value of noise intensity. 
Let us discuss these regimes in detail. 

A. Low noise intensity regime 

In this regime, instead of an exponential dependence 
on l/D, the curves shown in Fig. 2 display a power-law 
asymptotic behavior of the mean escape time 

T ^(a,D) = ^. (12) 

Further analysis of the data plotted in Fig. 2 allows us 
to determine the dependencies of the power law expo- 
nent 11(a) and coefficient C(a) on the Levy index a (see 
Fig. 3). The exponent n(a) is approximately unity up to 
a w 1.5. The dependence of C(a) for the first two po- 
tential types possesses a weak inflection at small a. Note 
that for fixed noise strength D, the dependence of T esc 
on the Levy index a is monotonic. 

Let us now turn to the probability density function 
(PDF) p(t) of the escape times. We use again our simu- 
lation scheme to obtain p(t) from the Langevin equation 
(11), but now for each fixed values of a and D we collect 
200,000 escape events, which are not averaged but pro- 
cessed with a simple routine, that constructs the PDF. 
The results for the three potentials from Fig.l are shown 
in Fig. 4 [34]. It is imminently clear that in all three 



cases and for all a (including the Gaussian case, a = 2) 
the exponential decay pattern 

p(t) = ±exp(-|) (13) 

is nicely obeyed, in agreement with the findings from 
previous simulations studies [23, 28, 29], as well as with 
the analytical results reported in Ref. [26] . 

Relation (13) allows us to extract the mean escape time 
T csc independently from our previous results in Fig. 2, in 
two different ways. Namely, we define 

T -=W) sTli (14) 

alternatively, we determine 

These two ways to determine T esc are intrinsically dif- 
ferent, T\ depending on the extrapolation of the fitted 
exponential behavior to t = 0, whereas T 2 is determined 
through the slope in the logarithmic versus linear plot. 
Tab. I shows the results for T\ and T 2 along with the val- 
ues of T esc obtained directly from the simulations results 
in Figs. 2. Indeed, the value of T 2 is consistently closer 
to T esc than T\. Overall, however, the agreement is very 
good (better than 1.5 per cent). 

B. High and intermediate noise intensity regimes 

We consider these two regimes for the case of the trun- 
cated harmonic potential as example, compare to the 
bottom graph in Fig. 2. The two other cases, bistable 
and metastable potentials, demonstrate similar behav- 
ior. When the noise intensity D increases beyond the 
low noise limit, in which the curves exhibit a linear de- 
pendence and are almost parallel, the various curves for 
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Figure 2: Mean escape times as function of noise intensity 
D for the bistable (top), metastable (middle), and truncated 
harmonic potentials (bottom). 
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Figure 3: Dependencies of (i(a) and C(a) for the bistable 
(top), metastable (middle) and truncated harmonic (bottom) 
potentials. See text. 



the escape time show a tendency to collapse. This is 
the intermediate noise intensity regime, which for the 
truncated harmonic potential lies approximately in the 
range —2 < —IgD < 1. Here, we observe that T csc 
demonstrates very weak dependence on the Levy index, 
in comparison with the low intensity regime. The ex- 



planation is the following. Contrary to the weak noise 
intensity regime, where the tail of the Levy noise distri- 
bution plays a decisive role in the escape process (less 
frequent but high amplitude noise spikes, that is, govern 
the barrier crossing), here the main role in the process, 
due to the higher intensities, passes on to the central part 
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Figure 4: Escape time probability density functions in depen- 
dence of time for the bistable (top) , metastable (middle) and 
truncated harmonic (bottom) potentials. In the log versus lin 
plot, the exponential dependence is obvious. 



of the distribution. Since all stable distributions, includ- 
ing the Gaussian, have central parts which are of simi- 
lar shape, the a-dependence of the escape time becomes 
weaker. Furthermore, in the intermediate region we ob- 
serve a non-monotonic dependence of T esc versus a, ex- 
amples for three values of D taken from the intermediate 
regime are shown in Fig. 5. This weak non-monotonicity 
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Table I: Mean escape times, obtained using three separate 
methods: T csc is determined directly from the simulations 
producing Fig. 2, while Ti and T2 are defined in Eqs. (14) 
and (15), respectively. 



most likely reflects minor differences in the shapes of the 
distribution of the noise £ a in the central part (relative 
to the much more pronounced differences in the tails of 
these PDFs) . This point will be investigated in more de- 
tail in a forthcoming paper. Another important feature 
of the intermediate regime is that the PDF of the escape 
times remains exponential. This is clearly seen in Fig. 6. 
In that figure, one can also recognize that for a = 0.1 and 
1.25, as well as a — 0.75 and 1.0 are almost coinciding. 

In Fig. 5 we also show the dependence of the mean 
escape time on the time increment St used in our simu- 
lations, showing a much higher sensitivity to the partic- 
ular value chosen for St than in the regime of low noise 
strength. This is intuitively clear, as now the process is 
terminated after few jumps only. As seen in Fig. 5, for 
increasingly small St, the shown curves converge. Note, 
however, that the weak non-monotonicity of the escape 
time on a is preserved for all St. 

When the noise intensity D is increased further, the 
intermediate regime turns over to the high noise intensity 
regime. Here, all the curves for the escape time T esc 
collapse into a single curve, that approaches the constant 
value T C sc = 10~ 2 , independent of D. Noting that in 
our simulation scheme the time step St = 10 -2 is used, 
we conclude that in the region of high noise intensity 
the particles reach the boundary in a single jump. We 
confirmed this conclusion by taking different values of the 
time step, observing that, whereas the picture in the low 
and intermediate regions changes negligibly, in the high 
intensity regime the collapsed curves again approach the 
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Figure 5: Dependence of mean escape time T csc on Levy index 
a, for various noise intensities D in the intermediate regime. 
For these results we used the truncated harmonic potential. 
The weak dependence on the time increment 8t of the simu- 
lation is discussed in the text. 



value of T C sc equal to a single time step. 




Figure 6: Escape time PDF for the intermediate noise 
strength, with D = 1. Note that the curves for a = 0.1 
and 1.25, as well as a = 0.75 and 1.0 are almost coincide, 
pointing at a weak inversion of the a-dependence. 



IV. ANALYTICAL RESULTS FOR THE 
CAUCHY CASE, a = 1 

A standard approach to the classical barrier crossing 
problem in the presence of Gaussian white noise is based 
on the stationary flux approximation assuming that the 
probability current across the barrier is constant. This 
is equivalent to requiring that the barrier is high in com- 
parison to thermal energy, or the noise intensity D low. 
The stationary flux approximation is widely used for the 
classical Kramers problem [1, 3, 6, 7, 35]. Here, we ex- 
tend this assumption to the case of white Levy noise, and 
obtain analytical results for the Cauchy case, a = 1, in a 
bistable potential. 

We start from the fractional Fokker-Planck equation 
(FFPE) 



9/ 
dt 



dx \ dx mj 



D 



d a f 
d\x\ a 



(16) 



for the density f(x,t) to find the diffusing particle at 
position x at time t [15]. Here, the potential U{x) is 
defined in Eq. (9a), and the fractional Riesz derivative 
d a /d\x\ a is understood via its Fourier transform, 



d\x\ 



:f(x,t)} = 



d\xj' 



■f(x,t)dx^~\k\ a f(k 1 t). 

(17) 



We note that the FFPE (16) is equivalent to the Langevin 
approach in Eqs. (4) to (7). 

After rescaling in the same way as outlined above, we 
arrive at the FFPE with dimensionless variables, 



W = d_fdU f 
dt dx V dx 



D 



d a f 
d\x\ a ' 



(18a) 
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where 



2 4 

f(*) = -y + V 



We can rewrite Eq. (18a) in terms of the flux j, 



df dj n 
— + — = 0. 

dt dx 



(18b) 



(19) 



We follow the stationary flux approximation of the Brow- 
nian theory, and also pursue the assumptions used to ob- 
tain the escape time in the Brownian limit [35] , and can 
be justified in a more stringent approach [3]. Namely 
we solve the stationary Eq. (16), under the assumption 
that the stationary distribution of particles differs signif- 
icantly from the equilibrium distribution: the "source" 
ensures that all particles are on the left-hand side of the 
barrier, causing the steady influx of particles toward the 
barrier; after crossing it, the particles disappear through 
the 'sink'. That is, following Ref. [35] we assume 



/" 

J — ( 



f(x)dx = 1, 



and 



j(x) = jo = T~ t 



(20a) 



(20b) 



where jo is the constant flux of particles across the bar- 
rier, and Tesc corresponds to the mean first passage time. 
Transforming to Fourier space, the FFPE becomes 

% = k m +k m- D w af w> (21) 



while Eq. (19) attains the form 



(22) 



Comparing Eqs. (21) and (22), we obtain for the Fourier 
transform of the flux the expression 

3(k) = + ^Dsigr ^ (k)\k^ 1 f. (23) 

Now, with the use of Eq. (20b) it follows from Eq. (23) 
that the stationary solution in the constant flux approx- 
imation is determined by the following equation [36]: 



dk 3 dk 



Dsign(k)f = 2mj S(k). (24) 



We solve Eq. (24) on the right and left semi-axes, and 
then match the solutions. The details are presented in 
Appendix B. This procedure yields 



/(*) 



jo2f&exp(z*fc) , fc>0, 
3o^exp(-zk), k < 



(25) 



where the constants a and b, and z are given in 
Eqs. (B21). The inverse Fourier transform is indeed a 
PDF, namely 



-ikx 



f(k) 



dk 
2^ 



J0^ j e k( iX +z) dk 

2ab 



Jo 



1 



2b [x + b) 2 + a 2 ' 



(26) 



The last step is the normalization conditions given by 
Eq. (20a). After integrating Eq. (26), 



2^(j+arctanQll= 1 



(27) 



we arrive at the final expression for the mean escape time, 

(28) 



* L 2 b 
T csc = -r-r 1 + - arctan - 
4ao \ tt \a 



For O « 1, we obtain from Eqs. (B26) and (B27) the 
asymptotic result for the escape time 



T csc ^-, £>« 1. 



(29) 



This result agrees with numerical simulations using the 
Levy noise £ Q as stochastic force, within the accuracy of 
about 12 per cent. 



V. CONCLUSIONS 

We consider symmetric Levy flights within the entire 
domain of Levy indices, < a < 2, in three generic types 
of the potential wells: bistable potential, metastable po- 
tential, and truncated harmonic potential. As the basis 
for the numerical analysis we use the Langevin equation 
with the Levy noise as source for the stochastic force, 
whereas for the analytical treatment we employ the space 
fractional Fokker-Planck equation. We obtain the follow- 
ing results. 

Firstly, we demonstrated by extensive numerical anal- 
ysis based on solving the Langevin equation the existence 
of three dynamic regimes in the barrier crossing behavior 
of a particle driven by Levy noise, namely: 

(i) The regime of low noise intensity D, displaying 
an algebraic dependence of the mean escape time on D 
through T csc ~ C(a)/D^ a K In this regime the nature of 
the Levy noise with its characteristic property to allow 
for large outliers, plays an essential role and dominates 
the escape process. We showed that the exponent /x(a) 
remains approximately constant, (i w 1 for < a < 2; 
close to a = 2, it displays a divergence towards the expo- 
nential dependence on 1/D at a = 2, the case of Gaus- 
sian noise. In this low noise intensity regime, we observe 
a monotonous increase of the escape time with increasing 
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a (keeping the noise intensity D fixed). The probability 
density of escape times decays exponentially 

(ii) The regime of intermediate noise intensity in which 
the escape time is determined by the central part of prob- 
ability distribution of the noise. This regime is charac- 
terized by the following features: The difference between 
the escape times for different Levy indices significantly 
decreases, since in the central part all stable distributions 
are very similar to each other, as well as to the Gaussian 
distribution. Additionally, there is a non-monotonic de- 
pendence of the escape time with increasing a. However, 
the probability densities of escape times still decay expo- 
nentially. 

(iii) The universal regime of high noise intensity where 
the particle escapes with the first step (or in very few 
steps) . The curves denoting the dependence of the escape 
time Tosc on noise intensity D collapse onto a single curve 
for all values of a (including the Gaussian limit). 

Secondly, for the particular Cauchy case, a — 1 we 
develop the kinetic theory of the escape over the barrier 
in a bistable potential. We start from the space frac- 
tional Fokkcr-Planck equation and use the assumption of 
a constant flux over the barrier. We find analytically the 
expression for the escape time, which at low noise inten- 
sities produces the result T csc w tt/D. This result agrees 
with numerical simulation within the accuracy of about 
12 per cent. 
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Then 



P{0<X<x) = P(0 < X < x, 7 > 0) 

=pfo<r^) v °-'<„>o N 



(W > x-"/( 1 - a )a(7), 7 > o) 



= P 
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The latter expression, according to Ref. [38] is an integral 
representation of Levy distribution. 

In the case 1 < a < 2 analogous steps pertain, but 
now for the quantity P (X > or, 7 > 0). When a = 1, 
the relation for X turns into X = tan 7, that, is known 
to be a random value with Cauchy PDF. One does not 
need to consider the case 7 < 0: it merely corresponds 
to negative X values. 

The numerical simulation of T esc as function of D is 
performed according to the following flowchart: 

1. A "particle" is placed at the potential minimum; 

2. A value of a is fixed; 



Appendix A: SIMULATIONS DESCRIPTION 

In this Appendix, we briefly outline our simulations 
strategy, in particular, how we simulate the white Levy 
noise £ Q (n). The corresponding generator is taken from 
[37]. We start by calculating the value 



X = 



sin (cry) / cos ([1 — a]j) 



(cos 7) 



W 



(l-a)/a 



(Al) 



where 7 is a random value uniformly distributed on the 
interval (— n/2, tt/2), W is an independent random vari- 
able with exponential distribution, and a is the Levy 
index, < a < 2. 

The proof that the calculated variable X possesses a 
Levy distribution function goes as follows. Firstly, let 
< a < 1. While 7 > 0, the expression for X may be 



3. The discrete-time Langevin equation (8) is iter- 
ated, until the particle reaches a definite coordi- 
nate, namely x = for the bistable potential, 
x = 10 for the metastable potential, and x = ±1 
for the truncated harmonic potential; 

4. The event time t = nAt is stored; 

5. Steps 3 and 4 are executed 10,000 times for each 
fixed value of D, and the average escape time cal- 
culated; 

A similar procedure is used to simulate the escape time 
PDF p(t), apart from taking the average. Explicitly: 

1 . The "particle" is placed at the potential minimum; 

2. Some value of D is fixed; 

3. Some value of a is taken; 



4. The discrete-time Langevin equation (8) is iterated, 
until the "particle" reaches a definite coordinate, 
namely x — for the bistable potential, x = 10 
for the metastable potential, and x = ±1 for the 
truncated harmonic potential; 

5. The event time t — nAt is stored; 

6. Step 4 is executed 200,000 times for each fixed value 
of a, but, this time, the obtained event times are 
not averaged, but handled with a simple routine, 
that builds the PDF; 

Each run was repeated and both results compared. 
Throughout, high reliability was observed. 



Appendix B: MATCHING PROCEDURE 

In this Appendix, we collect the necessary steps to 
identify the constants for the determination of the den- 
sity / in the barrier crossing for the Cauchy case. 

To this end, we solve Eq. (24) on the right and left 
semi-axes, and then match the solutions. For the two 
domains, we make the exponential ansatz 



where 



with 



f(k) = I Cl£Zk + ° 2eZ ' k > fc>0; 

\ C 3 e^ k + C^" k , fc<0, 



u + +v + + .u + -v + ^ 



3 D 

"+= 2 



3 D 

V += 2 



1 + 1/1 + 



27D 2 



1-1/1 + 



27D 2 



Moreover, we define 



+zV3 — 



with 



3 D 

U ~= 2 



■1 + 1/1 + 



27D 2 



v_ 



D 
~2 



-1-1/1 + 



27D 2 



= — ii 



+ ■ 



which implies that 

u+ + v + 



c 



-iV3 , 



<■ 



(Bl) 

(B2) 

(B3a) 
(B3b) 

(B4) 

(B5a) 
(B5b) 

(B6) 



To determine the unknown complex constants 
C\,C2,Cz,C/t we match the results at k = 0. Thus, 



firstly, we require f(k — 0) = f(k + 0), which, using 
Eq. (Bl) yields 



C\ + C*2 = C3 + C4. 



(B7) 



The second condition is obtained by integration of 
Eq. (24) over the small region [— e,e]: 

J dkf"' + J dkf'-Dj dksign(k)f = 2mj . (B8) 

In the limit £-*0we recover the condition 

/"(0+) - /"(0-) = 2nij , (B9) 
or, after inserting Eq. (Bl) into Eq. (B9), 

dz 2 + C 2 z* 2 - C 3 ( 2 - C 4 C* 2 = 2tt IJo . (B10) 
The third condition is that the PDF is a real function: 



Jo 



Hk f° rlk 

/(■'■)= / -z~ tkx h(k) + J ^e- ik *f 2 (k). (LSI 1 ) 
and 



2tt 



r°° dk f° dk 

f*(x) = ^ -e^/r(fc) + J ^e ik *f;(k). (B12) 

On the other hand, by changing k — ► — fc in Eq. (Bll), 
we find 

/° rlk f°° rlk 

J£e ikX M-Q + 1 ^h(-k). (B13) 



Since /(x) is a real function, we see from Eqs. (B12) and 
(B13): 

A*(fc) = H-k), f* 2 {k) = h(-k). (B14) 

Let us now insert Eq. (Bl) into Eqs. (B14). From the 
first equation we have 



C\e z k + C* 2 e zk = C 3 e-< k + C±e 



(B15) 



The second equation gives rise to the same result. Since 
z* = —( (see Eq. (B6)), we have 

C* = C 3 , C 2 * = C 4 . (B16) 

Moreover, from Eqs. (B16) and (B7) we see that 

(Ci + C 2 )* =^3 + ^4 = ^!+ C 2 ; (B17) 

thus, from Eqs. (B16) and (B17), we obtain the following 
relations: 

C\r — C 3 r, C 2 r — C4R, (B18a) 
C37 = — Cu, C41 — —C21, (B18b) 
C 2 i = —Cu, C41 = —C31, (B18c) 
C21 = C31 = -C u = -Cu- (B18d) 
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Therefore, we are actually dealing with 3 constants only, such that 
namely, C 1R , C 2 r and Cj: 

d - Cia + iC u C 2 = C 2R - id, (B19a) 
C 3 = da - id , d = d R + id ■ (B19b) 



Inserting Eqs. (B19) into Eq. (BIO) produces 

d R (z 2 -c 2 ) + dR ( z * 2 -e 2 ) 

+id (z 2 - z* 2 + C 2 - C* 2 ) = 2Trij . (B20) 
For convenience, we define 



Then 



u + + v + V3 



z = — a + «&, £ = a + «&; 



(B21) 



z^ - C = -4ia&, 
z* 2 - C* 2 - 4zafe, 
z 2 - z* 2 + c 2 - C* 2 = o. 



(B23a) 
(B23b) 
(B23c) 



With the use of Eq. (B23) we get from Eq. (B20): 



C 2 r — CiR 



2ab' 



(B24) 



Now, combine Eqs. (B19) and (B24), and insert into 
(B22) Eqs. (Bl): 



h(k) = (dn + id) e zk + (C 1R + H - id) e z ' k , k > 
f 2 (k) = (d R - id) e Ck + (c 1R + ^ + e C * fe , k < 0. 



(B25a) 
(B25b) 



r 



We are looking for the stationary (non-equilibrium) 
distribution; therefore we assume C\r = C'i = in order 
to retain in Eq. (B25) only those terms, that are propor- 
tional to the flux j - This leads us to the result (25). 

To calculate the D <C 1 limit of the characteristic 
escape time (29), we obtain the following limiting be- 
haviours. Thus, from Eqs. (B3) and (B21) we see that 



u+ 



v+ 



and 



a^D/2, 6w 1, 
'6 s 



arctan 



arctan 



D 



7T 

2' 



(B26a) 
(B26b) 

(B27a) 
(B27b) 



Appendix C: COMPARISON TO ANALYTIC 
RESULTS FROM REF. [26] 



We recall that in Ref. [26], the authors start from the 
formula 



(exp (ik£ t )) = cxp (-£(a)i|fc| Q ) . 



(CI) 



with < a < 2, where £t is the Levy stable process 
in the Levy-Khintchin representation. The constant £ is 
defined as 



f°° 1 

€(a) = 2 - 
Jo 



2r(l-a fira 
cos — 

o « v 2 



cosy 



dy = 



) (C2) 



for a strictly smaller than 2. In Ref. [26], the authors 
investigate the system driven by the Levy stable process 
eZt, where e corresponds to the noise intensity: 

(exp (ikeZt)) = exp (-£(a)te a \k\ a ) . (C3) 

Note that e is not completely equivalent to our noise in- 
tensity D. Indeed, in the main body of the paper we 
investigate the system as driven by the Levy stable pro- 
cess L t , with 



(exp (ikL t )) = cxp (-Dt\k\ a ) . 



(C4) 



In this Appendix, we show that our results are consis- 
tent with those obtained by Imkeller and Pavlyukevich 
from a different approach to the barrier crossing problem 
of a particle exposed to Levy stable noise [26]. 



Since L t = eZ t , from comparison of Eqs. (C4) and (C3), 
we see that 

D = £(a)e a . (C5) 
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Consider the exemplaric case of the truncated har- 
monic potential 



U(x) = 



x 2 /2, -1 < x < 1 
0, otherwise 



(C6) 



From our formula T csc ~ C{a)/D^ a \ assuming: that 
fi(a) w 1 (that we found to be nicely fulfilled for small 
and intermediate values of a, see Fig. 3), we find that 



C(a)/D. 



(C7) 



The result reported in Ref. [26] for this potential profile 
at small noise intensities behaves like 



a 

2e" 



(C8) 



see Eq. (24) from Rcf. [26], with their constants a and b 
chosen as unity, a = b = 1. 

Thus, if we require T csc = X, then by virtue of 
Eqs. (C7), (C8), (C5), and (C2), we find the relation 

C(a) = |£(a)=r(l-a)cos(^). (C9) 



r(i) = i, 


a — 


o, 


tt/2, 


a = 


1 


1/(2 -a)- 


-> oo, a — 


2" 



From this relation we find the asymptotic behavior of 
C(a), 



(CIO) 



that agrees nicely with our numerical results, see Fig. 3 
in the inset of the bottom graph corresponding to the 
truncated harmonic potential. 

Fig. 7 shows a direct comparison between the analyt- 
ical result from Ref. [26] and our numerical results, by 
enforcing the /z(a) = 1 equality. The agreement is ex- 
cellent. We also show that taking a variable fi(a) into 
account, only a slight deviation is obtained. 

For comparison, Fig. 8 shows the analogous results 
for the case of the bistable potential. While general 
agreement is quite good, for Levy index decreasing below 
a = 1 , an increasing deviation is observed. We note that 
according to the results from Ref. [26], the coefficient C 
is 7r at a = 1, in accordance with our result (29). 
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Figure 7: Comparison of the coefficients C(a) obtained by 
using numerical simulations with the analytical behaviour of 
£(a) (full line) derived in Ref. [26] , for the truncated harmonic 
potential. While the triangles denote our results for fi = 1 
fixed, the dashed line refers to our results with fJ,(a), where 
we permit an explicit dependence on a. 
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Figure 8: Comparison of the coefficients C(a) obtained by 
using numerical simulations with the analytical behaviour of 
£(q) (full line) derived in Ref. [26] , for the bistable potential. 
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